function impedance(CONSTS, plot_data)
    
    phi = 0;
    m_max = 20;

    I_full = func_I_full(phi, m_max, CONSTS);
    I_odn_pl = sum_I_odn_pl(phi, CONSTS);

    coef_sgs_to_si = 9e+11;
    Z = coef_sgs_to_si ./ I_full;
    Z_odn_pl = coef_sgs_to_si ./ I_odn_pl;

    if(plot_data)
        fprintf('Z = %e + %ei Om\n', real(Z), imag(Z));
        fprintf('Z_odn_pl = %e + %ei Om\n', real(Z_odn_pl), imag(Z_odn_pl));
    end

end